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ABSTRACT 


An  automatic  spectral-fitting  procedure  based  on  a  simple  least- 
squares  error  criterion  is  devised  for  determining  the  source  parameters  of 
earthquakes  from  surface  wave  data.  When  this  procedure  is  supplemented  by 
the  complex  cepstrum  technique  for  removing  the  spectral  modulations  due  to 
multipath  interference,  results  are  significantly  improved.  Rayleigh  and  Love 
wave  spectral  data  recorded  by  the  VLPE  network  and  the  NORSAR  long -period 
array  for  fifty-one  earthquakes  in  the  Eurasian  continent  are  processed  by  this 
procedure.  The  results  are  discussed  on  a  regional  basis.  It  is  found  that  the 
surface  wave  spectral  characteristics  may  vary  even  among  earthquakes  located 
in  essentially  the  same  epicentral  region.  However,  the  variations  are  not  ran¬ 
dom  but  follow  certain  distinct  patterns.  With  proper  data  coverage  the  focal 
depths  of  relatively  small  earthquakes  can  be  estimated  within  a  20-km  range 
by  employing  jointly  the  automatic  spectral-fitting  procedure  and  the  complex 
cepstrum  technique. 
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SECTION  1 
INTRODUCTION 

It  has  been  shown  that  the  amplitude  spectra  of  surface  waves 
observed  at  a  few  sites  can  be  used  to  give  reasonable  focal  depth  estimates 
for  earthquakes  whose  source  mechanisms  are  known  (Tsai  and  Aki,  1970a,  b; 
Mendiguren,  1971;  Canitez  and  Toksoz,  1972).  Some  people  have  voiced  ser¬ 
ious  doubts  about  the  validity  of  this  approach  (e.  g.  ,  McGarr,  1972),  Thus, 
a  thorough  evaluation  of  the  method  using  a  large  ensemble  of  earthquake 

data  has  been  undertaken  in  the  past  fifteen  months.  The  results  are  present¬ 
ed  in  this  report. 

Because  of  our  primary  concern  for  adapting  the  method  for 
small  magnitude  earthquakes,  we  have  to  start  with  the  assumption  that  (1)  no 
prior  knowledge  about  an  earthquake  source  mechanism  will  be  available,  and 
(2)  for  a  given  earthquake  reliable  surface  wave  spectral  data  will  be  available 
only  at  a  small  number  of  sites.  The  first  condition  requires  that,  besides 
focal  depth,  other  source  parameters  such  as  dip  angle,  slip  angle,  strike 
direction  and  the  seismic  moment  all  have  to  be  treated  as  unknown  variables. 
The  second  condition  rules  out  the  more  familiar  approach  of  using  radiation 
patterns  of  individual  frequencies  for  which  a  good  azimuthal  coverage  is 
essential.  As  an  alternative,  simultaneous  usage  of  Rayleigh  and  Love  wave 
amplitude  spectra  over  the  frequency  band  from  0.  02  to  0.  05  Hz  will  be  made. 
A  method  which  takes  into  account  the  conditions  mentioned  above  is  de¬ 
vised  by  using  a  least- squares  error  criterion.  A  detailed  description  of  the 
method  is  given  in  Section  11. 

Preliminary  tests  of  the  method  on  a  small  number  of  earth¬ 
quakes  were  carried  out  and  the  results,  as  previously  presented  in  the  first 


semi-annual  report  (Tsai,  1972a),  were  promising  enough  to  suggest  further 
evaluations  of  the  method.  Since  then,  we  have  embarked  on  an  effort  to 
apply  the  method  for  determining  the  source  parameters  and  foccl  depths  of 
a  large  number  of  earthquakes,  fifty-one  in  total,  located  in  the  Eurasian 
continent.  In  order  to  study  the  spectral  characteristics  of  surface -wave s 
on  a  regional  basis,  we  have  classified  these  events  into  nine  groups  accord¬ 
ing  to  their  epicentral  locations.  The  results  are  presented  in  Section  III. 

Remembering  that  our  method  employs  the  theoretical  spectra 
of  Love  and  Rayleigh  waves  calculated  from  a  layered  half  space  of  nondis¬ 
sipating  elastic  medium  to  compare  with  their  observed  counterparts,  it  be¬ 
comes  necessary  to  make  corrections  for  attenuation  due  to  anelasticity  and 
to  remove  the  spectral  modulations  due  to  interference  of  multipath  arrivals 
from  the  observed  spectra  before  comparisons  can  be  made.  For  the  attenua¬ 
tion  corrections  we  shall  use  smooth  curves  for  attenuation  coefficients.  Un¬ 
fortunately,  these  curves  are  not  well  determined  for  the  frequency  band  con¬ 
cerned,  namely  from  0.  02  to  0.  05  Hz,  especially  for  the  Eurasian  continent. 

We  have  obtained  an  attenuation  coefficient  curve  by  using  Rayleigh  waves  re¬ 
corded  at  Chiang  Mai,  Thailand  (CHG)  and  Toldeo,  Spain  (TLO)  from  six  earth¬ 
quakes  in  the  Moluccas -New  Guinea  area.  The  curve  is  subsequently  applied 
to  all  our  data.  Since  the  propagatioi  path  between  CHG  and  TLO  is  far  from 
typical  of  most  continental  paths  with  which  we  are  concerned,  we  have  made 
a  comparative  study  of  the  estimates  of  source  parameters  and  focal  depths 
using  different  curves  of  attenuation*  coefficients.  The  results  are  discussed 
in  Section  IV.  As  for  the  removal  of  the  spectral  modulations  due  to  inter¬ 
ference  of  multipath  arrivals,  it  has  been  demonstrated  in  our  earlier  report 
(Tsai,  1972a)  that  the  complex  cepstrum  technique  is  highly  suitable  for  this 
purpose.  We  have  applied  this  technique  to  some  of  our  data.  A  comparison 
is  made  between  the  outcomes  obtained  with  and  without  use  of  the  complex 
cepstrum  technique  and  is  presented  in  Section  IV. 
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The  method  is  designed  to  handle  surface  wave  spectral  data 
at  one  or  more  sites.  Limited  by  the  availability  of  data  at  the  time  when 
this  study  was  conducted,  all  except  two  of  the  fifty-one  earthquakes  have 
single-site  coverage.  The  remaining  two  have  data  available  at  three  sites. 
Thus,  it  is  possible  to  make  a  very  preliminary  investigation  on  outcomes 
of  multiple -site  coverage  as  compared  with  single-site  coverage.  The  re¬ 
sults  are  discussed  in  the  last  part  of  Section  IV. 

Finally,  Section  V  gives  a  summary  of  the  results  gained  to 
date.  We  shall  also  venture  to  draw  some  conclusions  from  what  is  accom¬ 
plished  with  regard  to  the  performance  of  the  method  employed  for  the  sur¬ 
face  wave  spectral  analysis. 
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SECTION  II 


AN  AUTOMATIC  SPECTRAL-FITTING  PROCEDURE  FOR 
DETERMINING  EARTHQUAKE  SOURCE  PARAMETERS 
AND  FOCAL  DEPTH  FROM  SURFACE  WAVES 

We  shall  attempt  to  determine  the  earthquake  source  para¬ 
meters  and  focal  depth  from  the  amplitude  spectra  of  surface  waves  obtained 
at  a  limited  number  of  recording  sites.  This  approach  is  particularly  suit¬ 
able  for  earthquakes  of  relatively  small  magnitudes  whose  source  parameters 
and  focal  depths  are  otherwise  difficult  to  determine.  In  terms  of  data  utility 
the  method  about  to  be  described  is  designed  to  take  advantages  of  the  in¬ 
creasing  availability  of  digital  seismograms  furnished  by  the  long-period  array 
stations,  such  as  NORSAR  and  ALP  A,  and  by  the  newly  installed  VLPE  Net¬ 
work.  Since  the  method  makes  use  of  the  theoretical  surface  wave  spectra  as 

its  core  element,  we  shall  start  with  a  brief  discussion  on  these  theoretical 
spectra. 

The  earthquake  source  is  represented  by  a  double  couple  which 
varies  as  a  step  function  in  time.  The  strength  of  the  double  couple  is  called 
seismic  moment  and  designated  as  M  .  The  Fourier  spectra  of  Love  and  Ray¬ 
leigh  waves  due  to  a  point  source  of  this  kind  buried  in  a  layered  half-space 
can  be  expressed  in  terms  of  normal  mode  solutions  for  respective  wave  types 
(Haskell,  1964;  Harkrider,  1964;  Ben-Menahem  and  Harkrider,  1964;  and 
Saito,  1967).  We  shall  follow  Saitoh  formulation  in  the  present  study.  The 
lengthy  expressions  for  computing  the  theoretical  surface  wave  spectra  due 
to  this  particular  case  of  source  representation  can  be  found  elsewhere  (Tsai 
and  Aki,  1970),  and  so  will  not  be  given  here.  It  is  sufficient  to  say  that  the 
spectral  values  cf  Love  and  Rayleigh  waves  of  a  given  frequency  depend  on  a 
total  of  six  parameters.  Two  of  the  six  parameters,  namely  the  dip  angle  & 
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and  the  slip  angle  A  ,  specify  the  orientation  of  the  double  couple.  The  third 
source  parameter  is  the  seismic  moment  M  which  represents  the  strength 
of  the  double  couple.  The  fourth  parameter  h  indicates  the  depth  at  which 
the  equivalent  point  source  of  the  earthquake  dislocation  is  located  in  the  lay¬ 
ered  half-space.  The  remaining  two  parameters  r  and  0  ,  specify  the  cyl¬ 
indrical  coordinates  of  the  observation  point  with  respect  to  the  epicenter.  6 
is  azimuth  angle  measured  counterclockwise  from  the  strike  direction  <p  of 

the  dislocation  plane.  Because  the  distance  parameter  r  appears  only  in  the 

-1/2  -nr 

geometrical  spreading  factor  as  r  and  in  the  attenuation  factor  as  e 

which  are  independent  of  the  source  characteristics,  it  can  be  conveniently 
handled  by  specifying  a  constant  value,  say  2  000  km,  for  it  in  computing  all 
the  theoretical  spectra.  In  the  meantime  observations  made  at  distances 
other  than  the  specified  one  would  have  to  be  equalized  by  making  proper  cor¬ 
rections  for  the  geometrical  spreading  and  attenuation.  Thus,  there  remain 
five  parameters  to  be  determined:  the  seismic  mement  M  ,  the  dip  angle  &  , 
the  slip  angle  A  ,  the  strike  direction  </>  ,  and  finally  the  focal  depth  h  .  One 
common  approach  to  this  problem  is  to  study  the  radiation  patterns  of  Love 
and  Rayleigh  waves  formed  by  spectral  amplitudes  of  a  single  frequency  ob¬ 
served  at  a  multitude  of  azimuths.  In  this  case  a  good  azimuthal  coverage  of 
recording  sites  is  essential.  Unfortunately,  for  earthquakes  of  relatively 
small  magnitudes  which  are  our  main  concern, it  will  be  extremely  difficult, 
if  not  impossible,  to  have  enough  azimuthal  coverage  to  yield  meaningful  rad¬ 
iation  patterns.  Hence,  this  approach  is  not  suitable  for  our  purpose.  As  an 
alternative,  we  shall  adopt  a  spectral  approach  which  makes  use  of  the  spectral 
shapes  of  both  Love  and  Rayleigh  waves  f  -rmed  by  individual  spectral  moduli 
of  a  number  of  successive  frequencies  within  a  iven  frequency  band,  say  from 
0.  02  to  0.  05  Hz.  In  this  case  we  can  take  maximum  advantage  of  those  good- 
quality  spectral  data  available  at  only  a  small  number  of  recording  sites. 

If  the  fault  plane  solution  is  available  for  an  earthquake  from 
observations  of  the  first  motions  of  body  waves,  it  is  possible  to  determine 
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three  of  the  remaining  five  parameters,  namely,  b  ,  A  and  <p  .  Like  using 
surface  wave  radiation  patterns,  this  approach  is  practical  only  if  there  is 
enough  coverage  of  recording  sites.  Our  primary  concern  with  relatively 
small  earthquakes  again  imposes  a  severe  limit  on  this  possibility.  Conse¬ 
quently,  we  shall  proceed  on  the  assumption  that  no  prior  knowledge  on  b  , 

A  and  <t>  will  be  available.  In  other  words,  all  the  five  parameters,  i.  e.  , 

£  ,  A  ,  <P  ,  h  and  M  will  be  treated  as  unknowns  in  our  approach  toward  in¬ 
terpreting  the  surface  wave  spectra.  A  detailed  description  of  our  method 
follows. 

The  amplitude  spectra  of  Rayleigh  and  Love  waves  from  an 
earthquake  are  obtained  by  using  the  fast  Fourier  transform  technique.  If 
necessary,  the  complex  cepstrum  technique  is  applied  to  remove  the  spectral 
modulations  caused  by  multipath  interference.  For  a  detailed  discussion  of 
this  technique  the  readers  are  referred  to  one  of  our  previous  reports  (Tsai, 
1972a).  A  sample  rate  of  one  sample  per  six  seconds  is  adopted  by  decimat¬ 
ing  the  original  time  traces  which  use  a  sample  rate  of  one  sample  per  two 
seconds.  Prior  to  transforming,  a  time  trace  of  any  length  is  filled  with  addi¬ 
tional  zeros  at  its  end  to  make  up  2048  time  points.  This  results  in  1025  spec¬ 
tral  points  at  frequencies  ranging  from  zero  Hz  to  the  folding  frequency  at  (1/12) 
Hz.  with  a  fixed  frequency  increment  of  (1/6x2048)  Hz.  For  the  present  study  we 
shall  focus  our  attention  on  the  spectral  behavior  in  the  frequency  band  between 
0.  02  and  0.  05  Hz  which  is  dictated  largely  by  the  instrumental  response  of  the 
long-period  seismographs  in  operation.  This  is  done  by  selecting  twelve  evenly- 
spaced  frequencies  ranging  from  (40/2048)  Hz  to  (106/2048)  Hz.  The  spectral 
values  at  these  frequencies  are  first  corrected  for  the  instrumental  response, 
the  geometrical  spreading  and  attenuation,  and  then  compared  with  their  theo¬ 
retical  counterparts  to  yield  estimates  on  the  five  earthquake  source  para¬ 
meters.  A  simple  least-squares  error  criterion  is  adopted  for  judging  the 
goodness  of  fit  between  the  observed  spectral  curves  of  Love  and  Rayleigh 
waves  defined  at  the  twelve  discrete  frequencies  and  their  theoretical  counter - 
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parts  computed  for  a  given  set  of  source  parameters.  The  curve -fitting  pro¬ 
cedure  is  described  below. 


R  L 

Let  X..  and  X..  denote  the  observed  Rayleigh  and  Love  wave 

,  .  1J  .  1J,  „  .th  .  .  .th  .  „  R 

spectral  amplitudes,  respectively,  for  1  frequency  and  at  j  site.  Y  and 

t  y 

Y  denote  the  corresponding  theoretical  spectra  for  a  seismic  moment  of  one 

dyne -cm.  The  sum  of  the  squares  of  residuals  £  is  defined  as: 


€  =  I  ?  \(XR.  -  MY?1)2  +  (X1?  -  MY^)2]  (1) 

j  i  L  ij  ij  ij  ij  J 

where  M  represents  the  seismic  moment  of  the  earthquake  source.  For  a 
R  L 

given  set  of  Y,.  and  Y..  which  are  functions  of  the  four  source  parameters 
ij  ij 

(b>  X ,  <t>  and  h),  the  error  count  is  minimized  by  finding  the  seismic  mo¬ 
ment  M  according  to  the  following  equation. 
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In  other  words,  for  a  given  combination  of  6  .  A ,  <p  and  h  ,  we  first  find  the 
optimum  seismic  moment  M  from  equation  (2).  This  is  followed  by  inserting 
the  resultant  M  value  in  equation  (1)  to  compute  the  minimum  error  count  for 
this  particular  set  of  source  parameters.  We  then  proceed  to  vary  the  four 
source  parameters,  b  ,  X ,  <P  and  h  ,  in  a  systematic  manner  so  as  to  cover  a 
wide  range  of  source  mechanisms  and  focal  depth.  In  the  mean  time  the  error 
count  €  is  used  as  a  convenient  criterion  for  judging  the  goodness  of  fit  be¬ 
tween  the  observed  and  the  theoretical  surface  wave  spectra  among  various 
sets  of  source  parameters.  In  the  end  we  shall  take  as  our  final  solution  the 
one  set  of  source  parameters  whose  theoretical  surface  wave  spectra  fit  most 
closely  with  their  observed  counterparts.  To  the  extent  that  the  prescribed 
layered  earth  model  provides  an  adequate  approximation  of  the  real  world, 
this  method  can  be  expected  to  provide  reasonable  estimates  of  the  earthquake 
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source  parameters  on  a  routine  basis  from  the  long -period  digital  seismo¬ 
grams.  Preliminary  tests  on  the  method  were  conducted  during  the  early 
part  of  this  project  and  the  results,  as  described  in  our  previous  report 
(Tsai,  1972a),  are  quite  encouraging.  Since  then,  we  have  undertaken  a 
more  thorough  evaluation  of  the  method  by  applying  it  to  the  surface  wave 
spectral  data  of  a  large  ensemble  of  fifty-one  earthquakes  taking  place  dur¬ 
ing  1971-1972  in  the  Eurasian  continent.  The  results  are  presented  in  Sec¬ 
tion  III.  Before  looking  into  the  outcome  of  individual  cases,  it  should  help 
it  this  point  to  provide  an  outline  of  the  procedure  through  which  the  final 
solution  is  found.  To  begin  with  the  range  and  the  increment  are  specified 
for  each  of  the  four  source  parameters  to  be  varied  as  follows: 


Parameters 

Range 

Increment 

Dip  Angle  ( J  ) 

60°-  90° 

10° 

Slip  Angle  (A) 

-90°^  90° 

30° 

Focal  depth  (h) 

0  -  125  km 

5  km 

Strike  direction  (0) 

0  180° 

10° 

The  strike  direction  0  is  measured  clockwise  from  the  North  direction.  In 
order  to  cover  every  one  of  the  specified  cases,  the  four  parameters  are  var¬ 
ied  in  sequence.  For  a  given  source  mechanism  defined  by  the  first  two  para¬ 
meters,  &  and  A  ,  the  focal  depth  h  is  increased  from  0  km  to  125  km  by  a 
fixed  increment  of  5  km.  For  each  focal  depth,  the  strike  direction  0  is  var¬ 
ied  from  0  to  180°  in  even  increments  of  10°.  At  each  combination  of  & ,  A  , 
h  and  0  ,  the  theoretical  Love  and  Rayleigh  wave  spectral  amplitudes  are 
computed  at  the  twelve  frequencies  mentioned  earlier  for  those  azimuths  in 
which  the  corresponding  spectral  data  are  available.  Both  the  theoretical  and 
the  observed  spectral  amplitudes  are  then  inserted  into  equations  1  and  2  to 
yield  an  estimate  of  €  and  M  ,  respectively.  At  each  increment  of  0,  if  the 
new  error  count  is  smaller  than  the  existing  error  count,  then  it  will  replace 
the  existing  one  for  next  comparison.  All  the  source  parameters  and  the 
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theoretical  spectra  associated  with  this  new  error  count  are  also  carried  over 
until  the  next  replacement  hy  an  even  smaller  error  count  takes  place.  In  this 
way,  we  arrive  at  the  optimum  strike  direction  0  lor  this  particular  source 
mechanism  at  a  particular  focal  depth.  This  same  procedure  is  applied  through- 

out  all  of  the  specified  focal  depths.  At  the  end  we  find  a  focal  depth  with  the 

"best-fit"  theoretical  surface  wave  <innrt„  r 

nace  wave  spectra.  For  every  combination  of  $  and 

we  keep  this  optimum  focal  depth  together  with  the  error  count  <  the 

strike  direction  0  ,  the  seismic  moment  M  and  the  theoretical  surflce  wave 

spectra.  There  are  a  total  of  twenty-eigh,  cases  of  source  mechanism  tested 

for  each  earthquake.  We  shall  look  into  the  distribution  of  the  corresponding 

twenty-eight  best  estimates  of  focal  depth.  By  so  doing  we  hope  to  establish 

a  probable  range  of  focal  depth  for  an  earthquake  regardless  of  its  source 

mechanism  because  knowledge  of  focal  depth  of  an  earthquake  can  be  very  use- 

ful  for  seismic  discrimination  purposes. 

Following  the  general  description  of  the  method  given  in  the  pre- 
ceding  paragraphs,  i,  is  appropriate  at  this  point  to  look  into  the  outcome  of 
actual  application  of  this  method  to  a  few  earthquakes.  For  this  purpose  we 
choose  six  earthquakes  that  occurred  in  the  Caucasus  region  during  1971-1972 
for  which  surface  wave  spectral  data  are  available.  The  PDE  information  for 
these  events  is  listed  in  Table  II- 1.  „  is  interesting  to  note  that  the  magni¬ 

tudes  of  these  events,  ranging  from  mb=  4.  0  to  mb=  5.  1,  are  relatively  small. 

Data  are  available  for  three  events  j 

s  with  code  names  starting  with  'LX'  at  KON 

a  VLPE  station  at  Konsberg,  Norway.  Data  are  also  available  a,  the  NORSAr’ 
long-period  array,  located  in  Norway,  for  the  other  three  events  with  their 
code  name  ending  with  'NIA  The  complex  cepstrum  technique  has  been  ap¬ 
plied  to  the  firs,  group  of  data  a,  KON,  bu,  not  for  the  second  group  a,  NORSAR. 
Following  is  an  examination  of  the  results  of  the  individual  events. 

Figure  II- 1  shows  the  results  for  Events  No.  17  and  18  which 
are  code  named  as  TRS/25 1 /22NE  and  CAU/283/09NL.  respectively.  Even, 
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PDE  INFORMATION  FOR  EARTHQUAKE -EVENTS  IN  CAUCASUS 


Amplitude  Spectral  Density  ( 


No.  17  is  located  near  the  Caucasus -Turkey  border  region.  The  spectra  of 
Rayleigh  and  Love  waves  recorded  at  NORSAR  are  shown,  respectively  in  the 
upper  left  and  upper  right  of  Figure  II- 1.  Both  of  the  observed  spectra,  as 
shown  in  dashed  curves,  are  contaminated  by  modulations  due  to  multipath 
interference.  The  solid  curves  represent  the  "best-fit"  theoretical  spectra 
which  result  in  a  minimum  error  count  €  .  The  corresponding  focal  depth  is 
5  km.  The  histogram  of  focal  depth  determinations  of  the  twenty-eight  cases 
of  source  mechanism  tested  for  this  earthquake  is  shown  in  the  lower  left 
of  Figure  II-l.  For  all  but  two  cases  of  source  mechanirm,  the  focal  depth 
is  set  at  0  or  5  km.  We  should  hasten  to  point  out  that  this  high  concentration 
of  focal  depth  estimates  in  a  narrow  range  does  not  necessarily  mean  a  high 
quality  of  focal  depth  determination.  This  can  be  seen  more  clearly  if  we  take 
a  look  at  Figure  II-2  in  which  the  error  count  is  plotted  against  the  correspond¬ 
ing  focal  depth  estimate  for  the  twenty-eight  cases  of  source  mechanism.  In 
this  figure  all  except  four  cases  have  nearly  identical  values  of  error  count 
which  simply  means  poor  resolution  for  favoring  a  particular  source  mechan¬ 
ism  based  on  the  error  count  criterion.  Such  phenomenon  is  to  be  expected 
for  highly  modulated  spectral  data.  Thus,  it  is  important  to  take  into  account 
the  quality  of  spectral  data  in  judging  the  quality  of  the  "best-fit"  solution. 

The  spectral  data  obtained  at  NORSAR  for  Event  No.  18  are 
shown  by  the  dashed  curves  in  the  middle  part  of  Figure  II-l,  with  Rayleigh 
waves  on  the  left  side  and  Love  waves  on  the  right.  The  corresponding  "best- 
fit"  theoretical  spectra  are  given  by  the  solid  curves.  The  histogram  of  focal 
depth  estimates  is  shown  in  the  lower  right  part  of  Figure  II-l.  All  focal  depth 
estimates  in  this  case  is  confined  in  the  upper  20  km  with  the  "best-fit"  focal 
depth  located  at  15  km.  The  error  count  is  plotted  against  focal  depth  in  Fig¬ 
ure  II- 3.  Judging  from  the  distribution  of  error  counts  and  the  goodness  of  fit 
between  the  observed  and  theoretical  spectra,  we  are  inclined  to  take  15  km  as 
a  reasonably  good  estimate  of  the  focal  depth. 
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Other  source  parameters  corresponding  to  the  best-fit  spectra 
for  the  preceding  two  earthquakes  are  listed  in  Table  II-2.  According  to  these 
parameters  Event  No.  17  is  an  oblique-slip  dislocation  with  its  strike  in  N140° 
E  whereas  Event  No.  18  is  a  normal  dislocation  on  a  vertical  plane  striking 
in  N30  E.  It  must  be  cautioned  that  these  dislocation  parameters  contain 
certain  ambiguities  due  to  the  symmetry  conditions  of  surface  wave  radiation 
from  a  double  couple  point  source.  These  ambiguities  can  not  be  resolved  by 
using  surface  wave  data  alone. 

The  observed  spectra  of  Rayleigh  (LR)  and  Love  (LQ)  waves  at 
KONfrom  Events  No.  19  and  20  are  shown  by  the  dashed  lines  in  Figure  II-4. 

As  indicated  in  Table  II- 1.  Event  No.  19  is  located  a  few  degrees  to  the  North¬ 
west  of  Event  No.  20.  These  spectra  have  been  processed  by  the  complex  cep- 
strum  technique  so  that  they  are  considerably  free  of  modulations.  The  spec¬ 
tral  behaviors  of  both  LR  and  LQ  waves  are  somewhat  similar.  The  best-fit 
theoretical  spectra  are  given  by  the  solid  curves  in  the  same  figure.  The 
agreement  between  the  observed  and  theoretical  spectra  is  generally  good  ex¬ 
cept  for  the  unsightly  presence  of  theoretical  LR  spectral  notches.  It  should 
be  emphasized  here  that  the  spectral  analysis  technique  for  determining  the 
source  parameters  of  an  earthquake  does  not  rely  solely  on  the  presence  of 
the  spectral  notches.  Instead,  the  technique  really  depends  more  on  the  gen¬ 
eral  shapes  of,  and  the  relative  levels  between,  LR  and  LQ  spectra.  The  source 
parameters  related  to  the  best-fit  theoretical  spectra  shown  in  Figure  II-4  for 
both  Events  No.  19  and  20,  as  listed  in  Table  II-2,  suggest  that  they  are  all 
associated  with  oblique- slip  dislocation;  one  being  40  km  deep,  and  the  other 
25  km  deep.  Realizing  that  the  method  being  used  for  the  present  study  can 
not  be  expected  to  yield  an  exact  estimate  of  focal  depth  for  an  earthquake,  it 
is  perhaps  more  instructive  to  set  a  probability  range  for  the  focal  depth  than 
to  pinpoint  it  at  a  particular  value.  For  this  purpose  the  histograms  of  twenty- 
eight  focal  depth  estimates  for  both  Events  No.  19  and  20  are  shown  on  Fig¬ 
ure  II-4.  For  Event  No.  19  the  estimated  focal  depth  is  confined  mostly 
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FIGURE  II-4 

LR  AND  LQ  SPECTRA  FOR  EVENTS  19  AND  20  (KON) 
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between  25  and  40  km.  A  study  of  the  profile  of  the  corresponding  twenty- 
eight  error  counts,  as  shown  in  Figure  II-5,  further  narrows  this  range  to 
between  30  and  40  km  in  which  many  relatively  small  error  counts  take 
place.  As  for  Event  No.  20,  the  focal  depth  estimates  spread  over  a  rather 
wide  range  from  5  to  35  km.  Again  a  glance  of  the  profile  of  error  counts, 
as  shown  in  Figure  II-6,  reveals  that  the  most  probable  focal  depth  range  for 
this  earthquake  should  be  between  15  and  35  km. 

The  LR  and  LQ  spectra  of  the  last  two  earthquakes  in  the  Cau¬ 
casus  area  listed  in  Table  II- 1  as  Events  No.  21  and  22  are  shown  in  the  up¬ 
per  part  of  Figure  II-7.  The  data  for  Event  No.  21  are  obtained  by  the  VLPE 
instruments  located  at  KON  and  are  processed  by  the  complex  cepstrum  tech¬ 
nic  ue.  The  data  for  Event  No.  22  are  obtained  at  NORSAR  and  are  not  pro¬ 
cessed  with  the  complex  cepstrum  technique.  A  comparison  of  the  spectra  be¬ 
tween  these  two  events  clearly  shows  the  spectral  improvement  brought  about 
by  the  complex  cepstrum  technique.  Once  again,  the  solid  curves  represent 
the  "best-fit"  theoretical  spectra  for  which  the  source  parameters  are  listed 
in  Table  II-2.  Both  events  are  shown  to  be  related  with  oblique -slip  disloca¬ 
tion  (  &  =  70°,  A  =  30°).  The  corresponding  focal  depth  is  determined  to  be  45 
km  and  55  km,  respectively,  for  Events  No.  21  and  22.  The  distributions  of 
focal  depth  estimates  for  the  twenty-eight  cases  of  source  mechanism  for  the 
two  events  are  shown  in  the  lower  part  of  Figure  II-7.  For  Event  No.  21  the 
estimated  focal  depth  ranges  from  25  to  50  km.  Further  examination  of  the 
associated  error  counts  shown  in  Figure  II-8  suggests  that  the  most  likely 
value  of  focal  depth  for  this  earthquake  lies  between  45  and  50  km.  With  re¬ 
gard  to  Event  No.  22,  we  find  a  small  range,  from  45  to  55  km,  of  focal  depth 
distribution  which  is  related  to  poor  quality  spectral  data.  The  error  counts 
in  this  case  tend  to  have  similar  values,  as  shown  in  Figure  II-9,  which  in 
turn  makes  it  difficult  to  reach  a  clearcut  solution  for  the  source  parameters. 
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FIGURE  II-7 

LR  AND  LQ  SPECTRA  FOR  EVENTS  21  (KON)  AND  22  (NORSAR) 


In  summary,  the  preceding  examples  of  six  earthquakes  in  the 
Caucasus  area  clearly  indicate  that  the  complex  cepstrum  technique  can  bring 
about  highly  desirable  improvement  on  the  quality  of  spectral  data  and  in  turn 
make  the  present  method  more  effective.  Judging  from  the  results  obtained 
by  applying  our  method  to  the  single -site  surface  wave  spectral  data  for  Events 
No,  19,  20,  and  21  to  which  the  complex  cepstrum  technique  was  applied,  it  is 
not  unreasonable  to  say  that  the  focal  depth  of  an  earthquake,  regardless  of  its 
source  mechanism,  can  be  determined  within  a  range  of  about  20  km.  It  should 
be  emphasized  that  f>  he  present  method  to  be  effective,  especially  for  re¬ 
latively  small  earthq_*„i<es  whose  surface  waves  are  expected  to  be  recorded 
with  sufficient  quality  only  at  a  very  small  number  of  sites,  simultaneous  use 
of  both  LiR  and  LQ  spectral  data  is  essential. 

We  have  applied  the  method  described  in  the  preceding  para¬ 
graphs  to  fifty-one  earthquakes  in  the  Eurasian  continent  including  the  six 
earthquakes  in  the  Caucasus  area  just  discussed  above.  The  spectral  data  of 
these  earthquakes  are  obtained  from  LP  seismograms  recorded  by  either  the 
Very  Long  Period  Experiment  (VLPE)  network  or  the  Norwegian  Seismic  Array 
(NORSAR).  The  VLPE  instruments  were  installed  and  operated  by  Lamont- 
Doherty  Geological  Observatory  of  Columbia  University.  A  copy  of  the  data 
tapes  was  made  available  to  the  Seismic  Data  Analysis  Center  (SDAC)  where 
routine  processing  and  analysis  was  performed  by  Texas  Instruments  Incor¬ 
porated.  A  detailed  description  of  the  computer  programs  written  for  this 
task  is  given  in  the  Documentation  of  Very  Long  Period  Experiment  Software 
Package  (Texas  Instruments  Incorporated,  1971a).  For  the  present  study,  only 
those  earthquakes  whose  surface  waves  have  sufficient  signal-to-noise  ratio  are 
selected  for  analysis.  The  VLPE  data  are  processed  by  the  complex  cepstrum 
technique  in  addition  to  a  bandpass  filtering  operation  for  the  pass  band  between 
0.  017  and  0.  070  Hz  to  remove  the  spectral  modulations  due  to  multipath  inter¬ 
ference.  Sixteen  out  of  the  fifty-one  earthquakes  presented  in  this  report  de¬ 
pend  on  the  VLPE  data  and  are  identified  by  event  codes  beginning  with  'LX'. 
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In  addition  each  event  code  contains  five  letters  in  the  middle  block  indicating 
the  geographical  location  of  the  earthquake  and  in  the  last  block  three  numer¬ 
ical  digits  for  bookkeeping  purposes  in  the  VLPE  routine  processing.  For 
example,  Event  No.  20  discussed  previously  is  given  an  event  code  LX/CAUCA/ 
606  to  indicate  that  the  VLPE  data  are  used  for  the  analysis  of  the  earthquake 
located  in  Caucasus  and  is  to  be  identified  as  Event  No.  606  in  the  VLPE  data 
series. 


Because  the  VLPE  network  was  still  in  experimental  stages 
during  the  period  when  this  study  was  undertaken,  the  amount  of  VLPE  data 
suitable  for  this  study  was  relatively  small.  As  an  alternative,  a  substantial 
quantity  of  the  NORSAR  data  was  taken  for  our  analysis.  This  was  conven¬ 
iently  accomplished  because  Texas  Instruments  Incorporated  was  also  res¬ 
ponsible  for  performing  routine  processing  and  analysis  of  the  NORSAR  data. 
Detailed  information  on  this  aspect  can  be  found  in  the  Documentation  of  the 
Array  Evaluation  Software  Package  (Texas  Instruments  Incorporated,  1971b). 
The  NORSAR  data  selected  for  use  in  this  study  were  the  be  am  steered  output 
of  surface  wave  trains  with  sufficiently  large  signal-to-noise  ratio.  The  com¬ 
plex  cepstrum  technique  was  not  applied  to  this  set  of  data  for  thirty-five 
earthquakes.  Each  event  code  for  the  NORSAR  data  is  set  to  begin  with  the 
first  three  letters  indicating  its  geographical  location  and  end  with  two  letters 
'NL'  identifying  itself  as  long -period  NORSAR  data.  In  between  there  are 
five  numerical  digits  showing  the  Julian  date  and  the  hour  in  which  the  event 
takes  place. 

The  results  of  our  analysis  of  these  fifty-one  earthquakes  will 
be  discussed  on  a  regional  basis.  In  total,  they  are  divided  into  nine  groups 
according  to  the  geographical  location  of  each  earthquake  epicenter.  For 
illustrative  purposes  the  results  of  six  earthquakes  in  Caucasus  already  have 
been  presented  as  a  group  in  the  preceding  paragraphs.  Results  obtained  in 
this  study  for  the  remaining  forty-five  earthquakes  which  are  divided  into 
eight  geographical  groups  are  described  in  the  following  section. 
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SECTION  III 

DATA  PRESENTATION  BY  REGION 


A.  OUTLINE  OF  PRESENTATION 

There  is  reason  to  believe  that  the  source  mechanism  for  earth 
quakes  may  be  related  to  regional  tectonic  structures.  Hence,  it  is  of  particu 
far  interest  to  examine  the  surface  wave  spectra  and  their  solutions  of  earth¬ 
quake  parameters  on  a  regional  basis.  For  this  purpose,  fifty-one  events  are 
grouped  into  nine  regions: 

Region  Event  Numbers 


1. 

Italy 

1-6 

. 

Greece 

7-10 

3. 

Turkey 

11-16 

4. 

Caucasus 

17-22 

5. 

Ir  an 

23-28 

6. 

Tadzhik 

29-33 

7. 

Sinkiang 

34-43 

8. 

Tibet 

44-48 

9. 

Northeastern  Asia 

49-51 

The  nine  regions  with  their  associated  events  are  identified  in  Figure  III-l. 

In  the  following  discussion,  the  nine  regions  (except  region  4  - 
Caucasus,  which  was  presented  in  Section  II)  are  examined  in  detail,  with  the 
observed  as  well  as  the  'best-fit'  spectra  presented.  Qualitative  judgements 
on  the  matching  of  the  'best-fit'  spectra  with  the  observed  spectra  will  be  left 
to  the  reader.  In  making  any  judgments,  one  should  bear  in  mind  the  proce¬ 
dure  used  in  determining  the  'best-fit*  spectra.  First,  the  least  squares  pro- 
cedure  involves  determining  the  minimum  residual  from  the  total  energy  in 


both  the  Rayleigh  and  Love  waves,  and  the  resulting  solution  is  not  determined 
by  fitting  the  Rayleigh  or  Love  wave  spectra  individually.  Secondly,  since  the 
earthquake  focal  depth  is  the  most  important  source  parameter  for  discrimina¬ 
tion  purposes,  our  discussions  oi  the  outcomes  from  the  least-squares  proce¬ 
dure  will  be  made  with  special  emphasis  on  the  focal  depth.  In  the  lower  part 
of  each  spectral  illustration,  the  distribution  of  focal  depths  obtained  from  the 
28  sets  of  dip  ar.d  slip  angles  is  grouped  and  plotted  as  the  number  of  appear¬ 
ances  versus  the  focal  depth.  The  shaded  group  in  this  plot  contains  the  mini¬ 
mum  residual  depth,  and  this  depth  along  with  the  other  source  parameters  is 
used  to  calculate  the  'best-fit'  spectra  shown  on  the  associated  spectral  plots. 

B.  EARTHQUAKES  IN  ITALY 

For  the  region  in  Northern  Italy,  six  events  have  been  examined. 
The  PDE  information  for  these  events  are  given  in  Table  III- 1 ,  with  the  theo¬ 
retically  determined  source  parameters  given  in  Table  III-2. 

From  the  shape  of  the  Fourier  amplitude  spectra  of  the  Rayleigh 
and  Love  waves,  the  Italian  events  can  be  placed  in  two  subgroups  (Events  1,2 
and  Events  3-6).  The  first  subgroup  (Events  1  and  2),  as  shown  in  Figure  III-2, 
is  characterized  by  a  prominent  trough  in  the  LR  spectra  between  0.  03  and  0.  04 
Hz.  Although  these  events  have  similar  spectral  shapes  for  both  the  Rayleigh 
and  Love  waves,  the  'best-fit'  theoretical  spectra  are  somewhat  different.  For 
Event  1,  the  node  at  0.  035  Hz  is  over-characterized  by  the  computed  Rayleigh- 
wave  spectrum.  A  vertical  strike-slip  fault  is  suggested  for  this  earthquake, 
with  the  focal  depth  in  the  range  from  5  to  35  km  having  a  minimum  residual  at 
25  km  (Table  III-2).  For  Event  2,  the  suggested  mechanism  is  also  a  vertical 
strike-slip  fault  with  the  focal  depth  in  the  range  from  5  to  40  km  and  having  a 
minimum  residual  at  30  km.  It  is  interesting  to  note  that  the  solution  for  both 
events  yields  a  strike -direction  of  N130°E. 
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Event  3  represents  a  multiple  station  fit  (Figure  III-3).  The 
'best-fit'  source  parameters  are  given  in  Table  III-2.  According  to  these 
source  parameters,  the  earthquake  is  related  to  oblique -slip  dislocation.  The 
focal  depth  falls  mostly  in  the  range  from  25  to  40  km  with  six  exceptions  at 
0  km.  More  discussion  on  this  event  will  be  made  in  the  next  section. 

The  amplitude  spectra  for  Events  4  and  6  are  shown  in  Figure 
III-4.  For  Events  4  and  6,  the  'best-fit'  dip  and  slip  angles  given  in  Table  III-2 
suggest  oblique-slip  dislocation  source  mechanisms.  The  focal  depth  estimates 
for  Event  4  predominantly  fall  in  the  range  of  30  to  45  km.  Both  Events  4  and  6 
have  nearly  the  same  residual-profile,  with  the  latter  having  a  focal  depth  in 
the  range  of  35  to  45  km.  The  strike  directions  for  these  two  events  are  1  1  0°E 
and  N10  E,  respectively,  which  probably  represent  conjugate  planes.  On  the 
basis  of  surface  wave  data  alone,  it  is  not  possible  to  make  a  choice. 

Spectra  for  Event  5  are  shown  in  Figure  III- 5.  The  residuals 
have  the  same  profile  as  those  for  Events  3,  4,  and  6,  with  the  focal  depth  in 
the  40-45  km  range. 

In  general,  we  can  group  Events  1  and  2  and  Events  3,  4,  5,  and 
6,  by  similar  spectral  shapes  and  similar  residual  profiles.  The  first  group 
has  a  strike-slip  source  mechanism  whereas  the  second  group  has  oblique  slip 
source  mechanisms.  Also,  the  first  group  seems  to  be  located  at  slightly 
shallower  depths. 

C.  EARTHQUAKES  IN  GREECE 

Four  events  have  been  studied  for  the  region  in  Greece  (Table 
HI- 3).  The  data  for  these  events  were  obtained  from  NORSAR.  The  best-fit 
solutions  are  given  in  Table  III-4. 

From  the  observed  spectral  shapes  of  the  Rayleigh  and  Love 
waves,  the  four  events  can  be  classified  into  two  subgroups;  Events  7  and  8  in 
Western  Greece  (Figure  III-6)  and  Events  9  and  10  in  Eastern  Greece  (Figure 
III- 7 ) . 
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FIGURE  IU-4 

LR  AND  LQ  SPECTRA  FOR  EVENTS  4  AND  6  (KON)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  9  AND  10  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL,  DEPTHS 


For  Event  7,  the  differences  of  the  residuals  are  small  and  one 
cannot  make  a  clear  estimate  of  the  dip  and  slip  angles.  The  focal  depth  could 
be  in  eif  r  the  range  from  5  to  10  km  or  from  35  to  40  km.  For  Event  8,  it 
would  seem  impossible  that  the  computet’  Love -wave  spectra  would  yield  a  node 
at  0.  035  Hz  for  the  theoretical  model.  The  best-fit  solution  yields  an.  oblique 
faulting,  with  a  focal  depth  in  the  20  to  45  km  range. 

For  Events  9  and  10,  the  residuals  are  distributed  for  5  to  10  km 
focal  depths  and  are  not  sufficiently  varied  to  make  a  clear  selection  of  the  earth¬ 
quake  parameters. 

None  of  the  Greek  events  were  processed  with  the  complex  cep- 
strum  technique,  so  the  solutions  should  be  taken  with  caution. 

D.  EARTHQUAKES  IN  TURKEY 

Six  events  have  been  studied  for  the  region  in  Turkey,  and  these 
data  sites  are  indicated  in  Table  III-5.  The  'best-fit'  earthquake  parameters 
are  given  in  Table  III-6, 

^or  Event  11,  the  'best-fit'  solution  suggests  an  oblique-slip 
source  mechanism.  The  focal  depth  is  found  in  the  range  of  5  to  10  km  (Figure 
HI- 8)-  The  observed  Rayleigh-wave  spectrum  for  Event  12  has  a  node  at  0.  047 
Hz  which  is  not  reproduced  by  its  'best-fit'  theoretical  counterpart.  The  mini¬ 
mum  residuals  are  about  equally  distributed  in  two  ranges  of  focal  depth,  5  to 
15  km  and  35  to  45  km.  Although  a  clearcut  choice  is  difficult,  the  shallower 

depth  range  is  slightly  preferred  because  of  the  spectral  similarity  between  this 
event  and  Event  11. 

Event  13  (Figure  III- 9)  has  its  minimum  residuals  narrowly  dis¬ 
tributed  in  the  range  of  45  to  55  km  focal  depths,  with  the  best  estimate  of  the 
focal  depth  being  50  km.  The  difference  in  both  Rayleigh  and  Love  wave  spectral 
shapes  between  this  event  and  Events  11,  12  is  remarkable. 
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LR  AND  LQ  SPECTRA  FOR  EVENT  13  (KON)  AND  THE  CORRESPONDING 
DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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Event  14  is  a  multiple  station  fit  (Figure  III- 1  0) ,  with  the  event 
information  obtained  from  LASA  bulletins.  The  minimum  residual  gives  a 
focal  depth  in  the  range  of  35  to  40  km,  The  best-fit  solution  suggests  an  ob¬ 
lique  slip  source  mechanism  with  the  strike  direction  either  N160°E  or  its 
orthogonal  direction. 

The  spectral  data  for  Events  15  and  16  as  shown  in  Figure  Hi¬ 
ll,  are  contaminated  by  modulations,  especially  for  Love  waves.  For  Event  15, 
the  residual  errors  are  narrowly  distributed  in  the  range  of  55  to  70  km,  but  a 
clear  selection  of  the  earthquake  parameters  is  difficult.  The  Love  waves  for 
Event  16  have  a  sharp  node  at  0.  038  Hz  which  probably  is  due  to  propagation 
effects.  The  residual  errors  are  distributed  equally  in  two  narrow  ranges,  5 
to  10  km  and  45  to  55  km,  making  a  clear  choice  of  focal  depth  difficult.  Ac¬ 
cording  to  the  best-fit  solutions,  both  Events  15  and  16  have  oblique-slip  source 
mechanisms.  In  view  of  the  poor  data  quality  these  solutions  should  be  taken 
with  caution. 

E.  EARTHQUAKES  IN  THE  CAUCASUS  (See  Section  II) 

F.  EARTHQUAKES  IN  IRAN 

Six  events  were  studied  for  the  region  in  Iran  (Table  III-7),  and 
the  source  parameters  obtained  for  these  events  are  given  in  Table  III-8.  In 
general,  the  Love  wave  amplitude  spectra  have  much  higher  levels  than  the 
Rayleigh  wave  amplitude  spectra  for  the  events  in  this  region. 

For  Event  23  (Figure  III -12),  the  nodes  are  not  consistent  between 
the  'best-fit'  theoretical  and  the  observed  Rayleigh-wave  spectra.  The  residuals 
from  the  twenty-eight  subsets  are  distributed  in  the  narrow  range  of  focal  depth 
from  45  to  50  km,  with  the  'best-fit'  solution  indicating  a  vertical  strike-slip 
source  mechanism.  The  residuals  for  Event  24  indicate  a  depth  range  of  70  to 
90  km. 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  15  AND  16  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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FIGURE  III- 12 


LR  AND  LQ  SPECTRA  FOR  EVENT  2  3  (NORSAR)  AND  EVENT  24  (KON) 
AND  THE  CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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Event  25  (Figure  III -13)  has  a  profile  of  residual-errors  which 
indicates  that  the  focal  depth  may  be  in  the  range  from  60  to  90  km  for  this 
earthquake.  The  best-fit  solution  for  the  strike -direction,  dip  angle,  and  slip 
angle  suggests  an  oblique-slip  source  mechanism. 

For  Event  26  (Figure  III- 14),  both  the  Rayleigh  and  Love  waves 
have  a  high  spectral  amplitude  level.  The  residual  errors  are  scattered  in  a 
wide  range  from  45  to  105  km  with  the  minimum  value  taking  place  at  80  km. 

A  strike-slip  source  mechanism  is  suggested  by  the  best-fit  solution.  The 
residual  errors  for  Event  27  are  distributed  in  the  narrow  range  of  60  to  70 
km  (Figure  III- 14).  According  to  the  best-fit  solution,  the  event  has  an  oblique* 
slip  source  mechanism.  Since  Events  26  and  27  are  located  at  the  same  place 
and  their  spectra  are  not  too  different,  the  difference  in  the  inferred  source 
mechanisms  is  probably  not  real  but  merely  reflects  poor  resolution. 

The  focal  depth  for  Event  28  (Figure  III- 15)  is  put  in  the  range 
of  55  to  70  km.  The  other  earthquake  parameters  suggest  an  oblique -slip  source 
mechanism  for  the  event. 

Examining  all  of  the  results  for  this  region,  it  can  be  seen  that, 
the  Love-wave  spectral  shapes  show  higher  amplitude  values  at  lower  frequen¬ 
cies  and  lower  amplitude  values  at  the  higher  frequencies  than  is  normally  ob¬ 
served  in  the  band  (0.  02  to  0.  05  Hz)  of  interest.  This  has  caused  our  auto¬ 
matic  spectral  fitting  procedure  to  place  relatively  large  focal  depths  for  these 
earthquakes. 

G.  EARTHQUAKES  IN  TADZHIK 

Five  events  recorded  at  NORSAR  have  been  studied  for  the  re¬ 
gion  in  Tadzhik  (Table  0-9),  and  the  estimated  parameters  for  these  events 
are  given  in  Table  III  - 1 0.  The  excitation  of  the  Love-waves  are  generally  low¬ 
er  than  those  for  the  Rayleigh  waves  with  the  exception  of  Event  33,  which  is 
located  to  the  northeast  of  the  other  four  events.  Generally,  the  spectral  shapes 
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FIGURE  III- 13 

LR  AND  LQ  SPECTRA  FOR  EVENT  25  (KON)  CHG)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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FIGURE  HI- 14 

LR  AND  LQ  SPECTRA  FOR  EVENTS  26  AND  27  (KON)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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for  both  the  Rayleigh  and  Love  waves  are  curved  up  from  the  lower  frequen¬ 
cies  to  the  higher  frequencies,  and  tend  to  give  relatively  shallow  focal 
depths. 

For  Eve.  t  29,  the  residual  errors  are  evenly  distributed  in  the 
focal  depth  range  of  30  to  50  km,  as  shown  in  Figure  III- 16,  and  other  source 
parameters  suggest  an  oblique -slip  source  mechanism.  The  focal  depth  for 
Event  30  is  located  in  the  10  to  30  km  range  with  the  greatest  possibility  at  20 
km.  The  other  source  parameters  of  the  best-fit  solution  suggest  a  dip- slip 
source  mechanism  for  this  earthquake. 

Because  of  the  variations  of  the  spectral  shapes  for  Event  31 
(Figure  III- 17),  a  good  estimation  of  the  source  parameters  is  not  possible. 

The  residual  errors  are  very  close  to  each  other  and  are  scattered  in  three 
spots,  with  the  focal  depth  at  0,  5  or  35  km. 

The  observed  amplitude  spectra  for  Event  32  (Figure  III -18)  are 
modulated,  with  the  modulation  more  severe  at  the  higher  frequencies.  The 
profile  of  residual  errors  places  the  focal  depth  in  two  ranges,  from  5  to  10  km 
or  from  30  to  40  km.  The  best-fit  solution  in  Table  III-10  picks  an  oblique - 
slip  source  mechanism  for  this  event.  For  Event  33,  the  focal  depth  lies  in 
the  range  of  30  to  45  km  with  the  best-fit  dip  and  slip  angles  defining  an  oblique- 
slip  source  mechanism.  The  best-fit  theoretical  spectra  follow  quite  closely 
the  general  shapes  of  their  observed  counterparts  for  Events  32  and  33.  Even 
so,  the  spectral  modulations  may  still  keep  the  resolution  for  the  best-fit  sol¬ 
ution  relatively  low. 

H.  EARTHQUAKES  IN  SINKIANG,  CHINA 

Ten  events  have  been  studied,  using  NORSAR  data,  for  the  region 
in  Sinkiang,  China  (Table  III -11).  The  earthquake  parameters  obtained  for  these 
events  are  shown  in  Table  III- 12.  In  general,  the  spectral  shapes  of  the  Rayleigh 
waves  are  characterized  by  nodes  in  the  0.  025  to  0.  035  Hz  band.  The  source 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  29  AND  30  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENT  31  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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FIGURE  HI- 18 

LR  AND  LQ  SPECTRA  FOR  EVENTS  32  AND  33  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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parameters  suggest  a  dip  angle  near  70°-  80°  and  the  strike  direction  is  ap¬ 
proximately  N90°E  for  this  region. 

For  example,  Event  34  {Figure  III- 19)  has  a  node  in  the  Rayleigh- 
wave  spectrum  in  the  0.  025  to  0.  035  Hz  band.  The  residual  errors  yield  a  focal 
depth  in  the  range  of  30  to  45  km,  with  the  high  frequency  Love  wave  energy  in¬ 
dicating  that  a  deeper  focus  is  very  unlikely.  Event  35  (Figure  III- 19)  yielded 
residual  errors  indicating  a  focal  depth  in  the  range  of  55  to  75  km  with  a  higher 
possibility  in  the  60  to  65  km  range.  The  best-fit  dip  and  slip  angles  define  a 
strike-slip  source  mechanism.  The  difference  in  the  Love  wave  spectral  shapes 
between  Events  34  and  35  is  remarkable. 

The  residual  errors  for  Event  36  (Figure  III-20)  indicate  a  focal 
depth  in  the  5  to  15  km  range,  with  the  dip  and  slip  angles  consistent  with  an 
oblique-slip  source  mechanism.  The  spectral  shape,  for  Event  37  are  similar 
to  those  of  Event  36.  The  focal  depth  is  estimated  to  be  in  5  to  15  km  range. 

The  best-fit  dip  and  slip  angles  suggest  a  dip-slip  source  mechanism.  Both 

Events  36  and  37  yield  a  consistent  strike  direction  of  N80°E.  The  difference 

in  the  best-fit  estimates  of  dip  and  sHp  angles  between  Events  36  and  37  is  questionable. 

Events  38  and  39  have  nodes  at  0.  035  Hz  in  the  Rayleigh  wave 
spectra  (f  igure  III-21).  The  residual  errors  yield  relatively  shallow  source 
parameters  in  Table  III -12,  Event  38  is  an  oblique -slip  fault  whereas  Event  39 
has  a  strike -slip  source  mechanism. 

Event  40  has  deep  nodes  in  both  the  Rayleigh  and  Love  wave  am¬ 
plitude  spectra  (Figure  III-22).  They  are  probably  due  to  propagational  effects 
rather  than  to  the  source  mechanism.  Residual  errors  place  the  focal  depth  in 

the  0-25  km  range.  For  Event  41  (Figure  III-22),  the  focal  depth  is  in  the  5  to  30 
km  range,  with  the  strike  direction  estimated  to  be  N90°E.  According  to  tne 
best-fit  dip  and  slip  angles  this  event  has  an  oblique-slip  source  mechanism. 

The  node  on  the  'best-fit'  theoretical  Rayleigh-wave  spectra  for 
Event  42  (Figure  III-23)  seems  to  be  over -characterized.  Residual  errors 
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FIGURE  III- 19 

LR  AND  LQ  SPECTRA  FOR  EVENTS  34  AND  35  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  36  AND  37  (NORSAR)  AND  THE 
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FIGURE  III-21 


LR  AND  LQ  SPECTRA  FOR  EVENTS  38  AND  39  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  40  AND  41  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 


No.  40,  LR,  NORSAR 


''N  * 


No.  4  0,  LQ,  NORSAR 

f' 


No.  4  L ,  LQ.  NORSAR 


No.  41  SIN/221/  0INL 


No.  40  SIN*064*04NL 


III-42 


Hz 


FIGURE  III-23 

LR  AND  LQ  SPECTRA  FOR  EVENTS  42  AND  43  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 


indicate  a  focal  depth  of  5  to  30  km  with  the  best  estimate  at  25  km.  The  Ray¬ 
leigh  wave  spectra  for  Event  43  also  has  a  deep  node,  which  seems  to  be  due 
to  path  effects.  Residual  errors  show  the  mechanism  as  an  oblique-slip  fault 
with  the  strike  direction,  N90°E,  consistent  with  the  other  events  in  this  re¬ 
gion.  The  distribution  of  focal  depth  estimates  spreads  over  a  relatively  wide 
range  from  0  to  4C  km.  The  best-fit  focal  depth  of  20  km  appears  to  be  rea¬ 
sonable  because  of  the  relatively  good  spectral  fitting  as  shown  in  Figure  111-2  3. 

Generally  speaking,  the  LQ  spectra  are  much  less  variable  than 
the  L.R  spectra  observed  at  NORSAR  for  the  earthquakes  in  this  region.  With 
the  exception  of  Event  35,  all  of  the  other  nine  earthquakes  in  this  region  are 
placed  at  relatively  shallow  focal  depths  in  the  upper  25  km. 

I.  EARTHQUAKES  IN  TIBET,  CHINA 

Five  events  for  the  region  in  Tibet  have  been  studied  by  using 
data  recorded  at  NORSAR  (Table  III -13).  The  best-fit  earthquake  parameters 
arc  given  in  Table  III- 14. 

For  Event  44,  (Figure  III-24),  the  observed  Love-wave  spectrum 
is  very  modulated.  The  profile  of  residuals  indicate  a  focal  depth  of  15  to  35  km, 
with  the  estimates  of  the  dip  and  slip  angles  suggesting  an  oblique-slip  faulting. 
The  spectral  fitting  for  Event  45  suggests  a  focal  depth  in  the  range  of  35  to  45 
km  and  an  oblique -slip  source  mechanism. 

The  observed  Rayleigh-wave  spectrum  for  Event  46  (Figure  III- 2 5 ) 
is  highly  modulated  and  at  a  lower  level  than  its  counterpart  for  Love  waves.  All 
residual  errors  are  close  to  each  other  with  a  maximum  to  minimum  ratio  of  1.  1  . 
The  focal  depth  may  be  in  the  range  of  5  to  40  km  with  the  greatest  possibility  at 
30  km,  and  the  dip  and  slip  angles  suggest  a  strike -slip  source  mechanism.  The 
Rayleigh-wave  is  very  modulated  for  Event  47  (Figure  III-25),  with  the  residual 
errors  confining  the  focal  depth  to  the  vicinity  of  45  km.  The  source  mechanism 
can  be  described  as  oblique  slip  faulting  with  the  strike  direction  N10°E  or  its 
orthogonal  solution  of  N100°E. 
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FIGURE  III -24 


LR  AND  LQ  SPECTRA  FOR  EVENTS  44  AND  45  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  46  AND  47  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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For  Event  48,  the  Rayleigh  wave  spectrum  is  also  highly  mod¬ 
ulated  (Figure  111-26),  with  the  residual  error  analysis  yielding  a  focal  depth 
of  2  5  to  40  km. 

In  general,  we  see  that  the  Rayleigh  wave  spectra  from  events 
in  Tibet  are  severely  modulated,  with  these  modulations  probably  due  to 
multipath  interference. 

J.  EARTHQUAKES  IN  NORTHEASTERN  ASIA 

The  last  three  events  studied  with  the  data  obtained  from  NOR- 
SAR  are  scattered  in  North  Eastern  Asia  (Table  III-  15).  The  source  mech¬ 
anism  solutions  are  given  in  Table  III-16.  All  the  observed  spectra  are  sev¬ 
erely  modulated,  as  shown  in  Figure  III-27  for  Event  49.  There  are  deep 
nodes  at  0.  03  Hz  in  both  the  Rayleigh  and  Love  wave  spectr  a,  and  hence  the 
solutions  are  not  considered  reliable. 

For  Events  50  and  51,  the  spectra  show  some  characteristics 
in  common  even  though  severely  modulated.  Thy  focal  depths  are  estimated 
to  lie  between  5  and  20  km,  as  shown  in  Figure  III -28. 
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LR  AND  LQ  SPECTRA  FOR  EVENT  49  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 
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LR  AND  LQ  SPECTRA  FOR  EVENTS  50  AND  51  (NORSAR)  AND  THE 
CORRESPONDING  DISTRIBUTION  OF  ESTIMATED  FOCAL  DEPTHS 


SECTION  IV 


EVALUATION  OF  THE  AUTOMATIC  INTERPRETATION 
PROCEDURE  FOR  SURFACE  WAVE  SPECTRA 

In  the  preceding  sections  we  have  made  a  factual  presentation 
of  the  results  obtained  from  applying  an  automatic  interpretation  procedure  to 
surface  wave  spectra  in  the  frequency  band  between  0.  02  and  0.  05  Hz  to  deter¬ 
mine  the  source  parameters,  including  focal  depth,  of  a  sizable  suite  of  fifty- 
one  Eurasian  earthquakes.  We  shall  now  turn  to  the  important  task  of  making 
a  thorough  evaluation  on  many  facets  of  the  procedure  which  may  have  vary¬ 
ing  degrees  of  influence  on  the  outcomes  derived  thereby.  Among  the  things 
to  be  covered  in  this  section  are  a  preliminary  survey  of  the  data  availability 
for  the  procedure,  the  effects  of  attenuation  corrections  and  complex  cepstrum 
technique  on  the  ol  -.ejved  spectra,  the  use  of  spectral  data  at  multiple  cites, 
and  a  summary  of  the  estimated  values  of  seismic  moment  for  the  fifty-one 
earthquakes  studied.  Let  us  now  proceed  to  discuss  these  topics  in  detail. 

A.  DATA  AVAILABILITY  FOR  THE  PRESENT  PROCEDURE 

The  automatic  interpretation  procedure  devised  in  this  study  is 
intended  for  making  maximum  use  of  the  surface  wave  data  of  relatively  small 
earthquakes  to  extract  information  useful  for  seismic  discrimination  purposes. 
Basically,  the  use  of  this  method  is  limited  only  by  the  availability  of  sufficient¬ 
ly  good  surface  wave  data.  Thus,  some  knowledge  of  the  magnitude  range  of  an 
earthquake  for  which  this  method  can  be  used  is  of  particular  interest.  Table 
IV- 1  gives  a  summary  on  the  mfa  magnitudes  and  the  epicentral  distances  at 
which  surface  wave  data  are  available  for  the  sixteen  earthquakes  studied  by 
using  the  VLPE  records.  Although  the  number  of  earthquakes  listed  in  the  table 
is  relatively  small,  it  appears  not  too  unrealistic  to  expect  that  for  many  Eur¬ 
asian  earthquakes  having  mb  magnitude  above  4.  0  a  sizable  portion  of  the  surface 
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wave  data  recorded  at  the  VLPE  sites  within  A  =  40°  will  be  usable  for  ana¬ 
lysis  by  our  method.  Hence,  the  amount  of  suitable  surface  wave  data  should 
increase  markedly  by  the  time  when  the  VLPE  network  is  in  full  operation. 

B.  EFFECTS  OF  THE  COMPLEX  CEPSTRUM  TECHNIQUE  ON  THE 
ESTIMATES  OF  EARTHQUAKE  SOURCE  PARAMETERS 

When  the  observed  and  the  best-fit  theoretical  spectra  for  both 
Rayleigh  and  Love  waves  for  the  six  earthquakes  in  the  Caucasus  area  were 
compared  in  Section  II,  we  have  pointed  out  that  application  of  the  complex  cep- 
strum  technique  on  the  VLPE  set  of  data  had  definitely  improved  the  quality 
of  the  observed  spectra  and  in  turn  resulted  in  a  marked  increase  in  the  re¬ 
solution  for  picking  the  correct  earthquake  source  parameters.  On  the  con¬ 
trary,  poorer  results  were  obtained  from  the  NORSAR  set  of  data  which  were 
not  processed  by  the  complex  cepstrum  technique.  This  phenomenon  was 
found  to  persist  throughout  the  data  for  the  other  forty-five  earthquakes  pre¬ 
sented  in  the  preceding  section.  Table  IV-2  compares  the  ratios  of  maximum 
to  minimum  error  counts  between  ten  earthquakes  in  the  VLPE  set  and  ten 
others  in  the  NORSAR  set;  all  are  randomly  chosen.  It  is  seen  that  the  ratios 
of  t  LPE  group  are  in  general  greater  than  the  NORSAR  group.  Since  the 
solutions  for  a  particular  earthquake  is  a  good  indicator  of  the  resolution  for 
deciding  on  the  best  solution,  we  can  say  that  the  source  mechanism  solutions 
for  earthquakes  from  the  VLPE  data  are  superior  to  the  NORSAR  data.  In 
summary,  it  is  appropriate  to  emphasize  the  improvement  brought  about  by 
the  complex  cepstrum  technique  both  on  the  excellence  of  spectral  fitting  and 
on  the  reliability  of  source  mechanism  solution  for  an  earthquake. 

C.  EFFECTS  OF  ATTENUATION  CORRECTIONS 

The  observed  amplitude  spectra  of  surface  waves  are  affected  by 
anelastic  attenuation.  Corrections  for  its  effects  have  to  be  made  before  the 
observed  spectra  can  be  used  for  determining  the  earthquake  source  parameters. 
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TABLE  JV-Z 

COMPARISON  OF  RESOLUTION  ERROR  BETWEEN  THE  VLPE 
DATA  {GEPSTR0M  APPLIED)  AND  NORSAR  DATA 
{NO  GEPSTRUM  APPLIED) 


VLPE  Data  (Ce’patrum  Applied) 


NORSAR  Data  {No  Gepeirum  Applied) 


Ratio  of  Residual 
(max/ min) 


Events 


LX  /  C  AU  C  A/  7  24 


LX/ITALY/ 665 
LX/SIR  AN/290 
LX/S1RAN/Z89 


LX/ SIR  AN/ 276 
LX/CITLY/134 


Ratio  of  Residual 

Events 

l  *  * 

3*28 

TUR/ 143/01 NL 

2.02 

TIB  /302/17NL 

1,32 

CAU/283/09NL 

1.74 

KRG/3Q1/13NL 

jh  -W.  U5S  v  | 

CAU/288/17NL 

3.94 

TRS/251  /22NL 

2.34 

TUR/120/04BE 

2.03 

TUR/221/04NL 

mA#0O6*O9NL 

2.39 

* 

1.66 

AIG*059*18NL 

Measurement  of  the  attenuation  coefficient  for  both  Rayleigh  and  JLove  waves 
in  the  frequency  band  between  0.  02  and  0.  05  Hz  have  been  made  by  a  number 
of  investigators.  Unfortunately,  these  measurements  were  made  either  as 
world-wide  averages  or  for  regions  outside  the  area  of  interest,  (i.  e.  ,  the 
Eurasian  continent)  and  hence  are  not  readily  applicable  for  the  present  study. 
To  improve  this  undesirable  condition,  we  have  made  an  attempt  to  measure 
the  attenuation  coefficient  of  surface  waves.  The  results  are  discussed  below. 

The  attenuation  coefficient  t}  (f)  for  a  given  frequency  f  is 

defined  as: 

-T)(f)(r  -  r  ) 

A2(f)  =  A^fje  (3) 

where 

Aj(f)  and  A„(f)  are  the  Fourier  amplitude  of  displacement 
for  frequency  f  at  Sites  1  and  2,  respectively; 
rj  and  r^  are  the  epicentral  distances  to  Sites  1  and  2,  respect¬ 
ively; 

T) (f )  is  the  attenuation  coefficient  for  frequency  f  . 

Equation  (  3  )  is  valid  only  if  the  seismic  source  and  Sites  1  and  2  are  aligned 
on  one  great-circle  path  since  in  this  case  the  effect  of  source  radiation  is 
identical  for  both  recording  sites. 

Six  earthquakes  with  magnitude  m  >  5.  0  in  the  Moluccas -New 

b 

Guinea  region  were  selected.  The  pertinent  information  about  these  events  as 
given  in  PDE  is  listed  in  Table  IV-3.  The  two  VLPE  sites  CHG  (Thailand)  and 
TLO  (Spain)  happen  to  be  aligned  approximately  along  the  same  great-circle 
paths.  Hence,  the  surface  wave  data  recorded  at  these  two  sites  for  the  six 
events  are  used  for  this  purpose.  The  Fourier  amplitudes  at  each  site  are 
corrected  for  instrumental  response  and  the  geometrical  spreading  before  they 
are  inserted  into  equation  (3)  for  computing  t)  (f ).  Because  the  data  quality  of 

horizontal  components  for  these  six  particular  events  is  not  sufficiently  good  for 
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the  present  purpose,  we  have  not  been  able  to  measure  the  attenuation  coeffic¬ 
ient  for  Love  waves.  Only  measurements  for  Rayleigh  waves  are  made  and 
reported  here.  The  attenuation  coefficient  11(f)  of  Rayleigh  waves  is  first 
computed  from  each  of  the  six  events  individually.  The  average  of  these  six 
values  is  taken  at  individual  frequencies  to  result  in  the  dashed  curve  shown 
in  Figure  IV- 1.  This  curve  is  further  smoothed  to  give  the  solid  curve  which 
has  been  used  to  correct  ail  the  observed  spectra,  Rayleigh  and  Love  waves 
alike,  presented  in  the  preceding  sections.  Our  measurements  suggest  that 
the  attenuation  coefficient  of  Rayleigh  waves  propagating  from  CHG  to  TLO 
increase  steadily  with  frequency  from  0.  02  to  0.  05  Hz.  This  behavior  is 
markedly  different  from  the  second  solid  curve  shown  in  Figure  IV-1  which 
was  measured  by  Tryggavasson  (1965).  Since  the  attenuation  correction  is 

(f)r 

implemented  by  multiplying  the  observed  spectral  value  by  a  factor  e 
at  any  given  frequency  f  ,  a  frequency  dependent  attenuation  coefficient  curve 
will  modify  the  original  spectral  shapes.  Apparently,  the  two  different  atten¬ 
uation  coefficient  curves  will  modify  the  same  original  spectral  shape  in  dif¬ 
ferent  ways.  Our  attenuation  coefficient  curve  will  result  in  an  increase  in  the 
spectral  slope  from  0.  02  to  0.  05  Hz  whereas  Tryggavasson's  curve  will  bring 
about  an  opposite  effect.  Thus,  the  source  parameters  estimated  from  the 
corrected  spectra  would  differ  somewhat  for  these  two  different  cases.  In 
order  to  see  how  much  difference  there  is  in  the  resultant  source  parameters, 
we  have  run  six  events  with  three  different  curves  for  attenuation  coefficient, 
namely,  our  present  curve,  Tryggavasson's  curve  and  a  zero  attenuation  coef¬ 
ficient  curve  for  the  whole  frequency  band.  This  last  curve  will  not  alter  the 
observed  spectral  shapes.  The  resultant  source  parameters  for  these  six 
earthquakes  are  listed  in  Table  1V-4.  In  the  table  ea"h  event  is  given  three 
different  solutions  A,  B,  and  C  which  are  obtained  according  to  the  three  cases 
of  attenuation  correction  just  mentioned  above.  The  results  suggest  that  Solu¬ 
tion  A  which  is  based  on  our  own  attenuation  coefficient  curve  has  a  tendency 
to  yield  a  shallower  focal  depth  than  Solution  B  which  uses  Tryggavasson's 
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Fortunately,  the  difference  stays  well  within  15  km  for  most  cases. 

In  other  words,  the  present  method  for  determining  focal  depths  is  not  criti¬ 
cally  dependent  on  precise  attenuation  corrections  so  long  as  the  attenuation 
coefficient  can  be  approximated  by  a  smooth  curve.  In  case  of  an  incorrect 
curve  for  attenuation  coefficient,  the  resultant  estimates  of  focal  depth  will  be 
biased  all  in  the  same  pattern.  This  will  not  alter  substantially  the  relative 
focal  depths  among  a  suite  of  earthquakes  which  are  of  particular  interest  to 
us.  Thus,  as  far  as  focal  depth  is  concerned,  lack  of  precise  information  for 
making  attenuation  corrections  on  the  observed  spectra  is  not  so  serious  as 
one  may  suspect.  It  must  be  cautioned  that  similar  remarks  can  not  be  extend¬ 
ed  to  other  source  parameters  such  as  the  strike  direction,  the  dip  and  slip 
angles  and  the  seismic  moment  because  substantial  differences  do  take  place 
on  the  estimates  of  these  parameters  from  different  solutions  as  can  be  seen 
from  Table  IV-4.  For  this  reason  it  is  strongly  hoped  that  more  definite 
measurements  on  the  attenuation  coefficients  of  both  Rayleigh  and  Love  waves 
can  be  made  on  a  regional  basis  for  regions  of  interest  in  the  near  future. 

D.  MULTIPLE-SITE  SPECTRAL  FITTING 

The  automatic  spectral  fitting  procedure  is  intended  for  handling 
surface  wave  data  at  multiple  sites.  Unfortunately,  the  applications  of  this  pro¬ 
cedure  to  the  earthquakes  described  in  the  preceding  two  sections  have  been 
made  mostly  on  single-site  data.  In  order  to  have  some  knowledge  on  the  per¬ 
formance  of  the  present  procedure  for  fitting  the  spectral  data  at  multiple  sites, 
we  have  made  a  study  of  two  earthquakes  for  which  data  are  available  at  more 
than  one  site.  The  two  earthquakes  are  Event  No.  3  (LX/ CITLY/ 1 34 )  in  Italy 

and  Event  No.  21  (LX/ CAUC A/724)  in  Caucasus.  The  results  are  presented 
below. 

Figure  IV-2  shows  the  observed  spectra  in  dashed  curves  for 
Event  LX/ CITLY/ 1  34.  Both  LR  and  LQ  spectra  are  available  at  KON  and  OGD 
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whereas  only  LQ  spectrum  is  usable  at  TLO.  All  these  three  sites  belong  to 
the  VLPE  Network.  The  distances  from,  and  the  azimuths  at  the  epicenter 
are  15.  5°  and  -7.4°  to  KON,  61.  6°  and  -59.  0°  to  OGD,  13.  7°  and  -102.  8°  to 
TLO,  respectively.  The  data  from  OGD  clearly  suggest  that  both  the  LR  and 
LQ  spectral  shapes  deteriorate  with  increasing  epicentral  distance.,  The  'best- 
fit  theoretical  spectra  are  superimposed  on  the  observed  spectra  in  Figure 
for  four  different  solutions  designated  as  A,  B,  C,  and  D,  respectively. 
Solution  A  is  obtained  by  using  the  observed  spectra  at  all  three  sites.  The 
fit  between  the  observed  and  the  theoretical  spectra  is  quite  acceptable  for  all 
three  sites.  The  corresponding  source  parameters  are  given  in  the  second  line 
of  Table  IV-5.  According  to  Solution  A,  this  earthquake  is  located  at  a  depth  of 

30  km  and  is  related  to  a  oblique  slip  dislocation  on  a  vertical  plane.  Solution  B 
is  obtained  by  using  the  observed  LR  and  LQ  spectra  at  KON  and  OGD.  The 
theoretical  spectra  of  this  solution  give  a  reasonable  fit  with  the  observed 
counterparts  at  OGD  but  are  too  low  to  explain  the  observed  spectra  at  KON. 

The  corresponding  source  parameters  are  given  in  the  third  line  of  Table  IV-5. 
The  focal  depth  is  estimated  at  25  km  in  this  case.  When  compared  to  Solution 
A,  Solution  B  appears  to  be  less  satisfactory  in  terms  of  the  goodness  of  spec¬ 
tral  fitting.  Solutions  C  and  D  are  obtained  by  using  the  single -site  spectral 
data  at  KON  and  OGD,  respectively.  The  theoretical  spectra  for  the  other  site 
are  then  'predicted'  according  to  each  single-site  solution.  For  Solution  C  the 
predicted  spectra  for  OGD  do  not  agree  with  their  observed  counterparts,  es¬ 
pecially,  the  LR  spectra.  Similarly,  for  Solution  D  the  predicted  spectra  for 
KON  are  much  too  high  to  fit  both  the  observed  LR  and  LQ  spectra.  The  cor¬ 
responding  source  parameters  for  these  two  solutions  are  given  in  the  last  two 
lines  of  Table  IV-5.  The  focal  depths  are  given  as  45  and  35  km,  respectively. 
In  this  table  we  find  that  the  strike  direction,  the  dip  and  slip  angles  as  well  as 
the  focal  depth  all  change  from  solution  to  solution  with  varying  degrees.  It  is 
difficult  to  make  a  definite  judgement  on  whether  the  multiple-site  coverage 
yields  more  reliable  estimates  on  the  source  parameters  than  the  single-site 
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coverage  using  the  good  quality  data  at  KON  alone,  because  the  quality  of  uata 
are  not  uniform  for  the  three  sites  considered  in  this  case.  Inclusion  of  poor 
quality  data  will  certainly  not  bring  about  improvement  on  the  quality  of  the 
solution  for  source  parameters.  Therefore,  data  quality  at  individual  sites 
should  not  be  overlooked  in  judging  the  multiple -site  solution.  This  point  is 
brought  out  clearly  by  the  next  example. 

Figures  IV-3  and  IV-4  show  the  observed  spectra  for  Event 
LX/CAUCA/724  at  the  three  VLPE  sites  FBK,  KON  and  OGD.  Only  the  LR 
spectrum  is  available  at  FBK  in  this  instance.  The  distances  from,  and  the 
azimuths  at  the  epicenter  are  30.  3°  and  -39.  2°  to  KON,  73.  9°  and  3.  8°  to  FBK, 
83.  2  and  -39.  4  to  OGD,  The  two  sites  KON  and  OGD  in  this  case  happen  to 
be  located  almost  on  the  same  great  circle  path.  Therefore,  the  spectra  ob¬ 
served  at  these  two  sites  should  be  very  similar  to  each  other,  if  they  are  not 
subject  to  grossly  different  effects  of  propagation.  A  comparison  of  the  obser- 
ed  LR  and  LQ  spectra  between  KON  and  OGD  clearly  reveals  that  the  spectra 
at  OGD  markedly  differ  from  their  counterparts  at  KON,  especially  for  the  LR 
waves.  This  shows  that  the  greater  the  epicentral  distance  becomes,  the  poor¬ 
er  the  spectral  quality  will  be.  Thus,  in  this  case,  addition  of  the  data  at  OGD 
to  the  automatic  spectral  fitting  procedure  would  actually  downgrade  the  out¬ 
comes.  Addition  of  the  poor  quality  data  at  FBK  will  have  the  same  effect  too. 
Now  let  us  look  at  the  'best-fit'  theoretical  spectra  and  the  related  source  para¬ 
meters.  Again  four  solutions  are  sought  and  shown  as  A,  B,  C,  and  D  in  Fig¬ 
ures  IV-3  and  IV-4,  respectively.  The  related  source  parameters  for  these 
solutions  are  given  in  Table  IV-6.  The  theoretical  spectra  corresponding  to 
Solution  A  which  is  obtained  by  using  the  data  at  all  three  sites  fail  to  give  a 
satisfactory  fit  to  their  counterparts,  especially,  for  the  LQ  spectra  at  both 
KON  and  OGD.  Solution  B  which  is  derived  by  using  data  at  KON  and  OGD  ap¬ 
pears  to  yield  somewhat  better  results.  Still,  it  fails  to  produce  an  acceptable 
spectral  fit  for  KON.  Solution  C  is  obtained  by  using  the  data  at  KON,  the  spec¬ 
tral  fit  is  excellent  at  this  site.  The  predicted  spectra  from  this  solution  again 
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fall  short  of  giving  a  satisfactory  spectral  fit  at  OGD.  As  for  Solution  D,  the 
outcome  as  even  less  desirable  because  it  fails  to  reproduce  large  enough  LQ 
spectra  at  both  KON  and  OGD.  The  focal  depth  estimates  in  this  case  are 
highly  variable  among  the  four  solutions.  In  summary,  the  addition  of  more 
data  with  poor  quality  in  our  automatic  spectral  fitting  procedure  definitely 
would  have  negative  effects  on  the  solution  for  source  parameters.  Thus,  one 
should  be  careful  during  the  selection  of  spectral  data  to  be  used  for  this  pur¬ 
pose.  Conversely,  we  can  expect  that  addition  of  more  data  with  good  quality 
would  improve  the  outcomes  of  the  present  procedure.  The  extent  of  this  im¬ 
provement  still  remains  to  be  evaluated  when  such  data  become  available. 

E.  ESTIMATION  OF  SEISMIC  MOMENT 

As  pointed  out  earlier  in  Section  II,  the  seismic  moment  is  used 
as  a  parameter  for  minimizing  the  error  in  our  spectral  fitting  procedure.  Thus, 
the  seismic  moment  is  determined  along  with  other  source  parameters  in  the 
process.  Figure  IV-5  shows  the  seismic  moment  corresponding  to  the  "best- 
fit''  theoretical  spectra  that  is  plotted  against  the  PDE  mfa  magnitude  for  each 
of  the  fifty-one  earthquakes  processed  by  the  present  method.  The  solid  line  in  the 
figure  represents  the  theoretical  curve  based  on  the  revised  -model  for  earth¬ 
quake  source  spectra  previously  proposed  in  our  second  semi-annual  tecnnlcal 
report  (Tsai,  1972b).  The  curve  was  calibrated  then  with  observations  for  earth¬ 
quakes  in  the  magnitude  range  higher  than  most  of  the  earthquakes  shown  in  the 
present  figure.  It  is  remarkable  that  the  theoretical  curve  agrees  so  closely  with 
the  data  obtained  in  the  present  study.  If  we  convert  mfa  to  M  magnitude,  this 
implies  that  the  seismic  moment  is  directly  proportional  to  the  M  magnitude  for 

M  smaller  than  5.  0. 
s 
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SECTION  V 


CONCLUSIONS  AND  DISCUSSIONS 

To  summarize  the  results  obtained  from  application  of  the  auto¬ 
matic  spectral-fitting  procedure  to  the  data  presented  in  the  preceding  sections 
the  following  remarks  are  in  order: 

•  The  use  of  a  relatively  large  number  of  unmodified  spectral 

data  in  this  study  has  been  intended  for  the  purpose  of  gaining 
first-hand  knowledge  on  the  spectral  behavior  of  surface  waves 
carried  by  the  original  recordings  for  earthquakes  in  regions 
of  interest.  It  is  found  that  a  sizeable  proportion  of  the  unmod¬ 
ified  spectral  data  are  severely  contaminated  by  modulations 
related  to  multipath  interference.  In  the  meantime,  the  oc¬ 
currence  of  this  phenomenon  is  found  to  be  highly  erratic.  In 
some  cases  the  spectral  modulations  take  place  on  both  Rayleigh 
and  Love  waves,  whereas  in  other  cases,  they  are  present  on 
only  one  of  these  two  types  of  surface  waves.  There  is  indica¬ 
tion  that  the  spectral  modulations  become  more  serious  with  in¬ 
creasing  epicentral  distance.  This  last  phenomenon  is  not  sur¬ 
prising  because  the  longer  the  propagation  path  becomes,  the 
more  lateral  heterogeneities  the  surface  wave  will  encounter. 
This,  in  turn,  will  result  in  more  severe  spectral  modulations. 
For  our  purpose  the  presence  of  such  undersirable  spectral 
modulations  will  degrade  the  outcomes  from  applications  of  the 
present  procedure  to  the  data.  Fortunately,  such  spectral  mod¬ 
ulations  can  be  removed  by  using  the  complex  cepstrum  technique. 
This  improvement  brought  about  by  the  technique  is  clearly  de¬ 
monstrated  in  the  sixteen  earthquakes  for  which  the  VLPE  data 
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are  processed  with  the  complex  cepstrum  technique.  Thus,  a 
combination  of  this  technique  with  the  automatic  spectral  fitting 
procedure  will  make  it  possible  for  us  to  use  the  surface  wave 
data  for  estimating  an  earthquake's  source  parameters  which 
may  not  be  available  otherwise. 

Our  primary  objective  has  been  centered  on  developing  a  pra¬ 
ctical  procedure  for  exploiting  the  long-period  surface  wave 
data  to  provide  us  useful  information  for  seismic  discrimination 
purposes  in  terms  of  not  only  the  M  magnitude  but  also  some 
other  source  parameters.  In  presenting  the  results  of  this  study 
we  have  singled  out  one  such  parameter,  the  focal  depth,  for 
special  emphasis.  The  rationale  behind  this  approach  is  the 
following:  For  a  given  set  of  surface  wave  data  if  we  can  esti¬ 
mate  the  focal  depth  within  a  relatively  narrow  range,  then  we 
shall  be  able  to  make  a  judgement  on  whether  the  seismic  event 
is  natural  or  man-made  based  on  its  location  in  depth.  Unfor  - 
tunately,  based  on  the  present  data  it  is  not  yet  possible  to  make 
a  definite  estimate  of  the  depth  range  attainable  by  the  present 
procedure.  It  should  be  emphasized  once  again  that  the  auto¬ 
matic  spectral  fitting  procedure  is  basically  limited  only  by  the 
availability  of  sufficiently  good  data.  There  is  reason  to  expect 
that  when  proper  multiple  -  site  data  are  available,  the  focal  depth 
of  an  earthquake  can  be  estimated  within  a  range  of  twenty  kilo¬ 
meters.  Further  investigations  on  this  aspect  are  needed. 

Since  our  present  procedure  depends  on  the  theoretical  surface 
wave  spectra  computed  for  a  layered,  nondissipating  earth  model, 
detailed  knowledge  on  the  elastic  and  anelastic  properties  of  the 
crust  and  the  upper  mantle  in  the  regions  of  interest  will  improve 
the  performance  of  the  automatic  spectral  fitting  procedure. 
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Preliminary  measurements  on  the  attenuation  coefficient  of 
Rayleigh  waves  have  been  made  for  the  path  between  CHG  and 
TLO.  Data  for  the  central  part  of  the  Eurasian  continent  re¬ 
mains  to  be  desired. 


SECTION  VI 
REFERENCES 


Ben  -Menahem,  A.,  and  D.  G.  Harkrider,  1964.  Radiation  Patterns  of  Seismic 
Suiface  Waves  from  Buried  Dipolar  Point  Sources  in  a  Flat  Stratified 
Earth,  J.  Geophys.  Res.,  6 9,  2605-2620. 

Camtez,  N.  ,  and  M.  N.  Toksoz,  ,972,  Static  and  Dynamic  Study  of  Earthquake 
Source  Mechanism,  J.  Geophys.  Res.,  77,  2583-2594 

Harkrider,  D.  G.  ,  1964,  Surface  Waves  in  Multilayered  Elastic  Media,  1, 

Rayleigh  and  Love  Waves  from  Buried  Sources  in  a  Multilayered  Elastic 
Halfspace,  Bull.  Seismol.  Soc.  America,  54,  627-679. 

Haskell,  N.  A.,  1964,  Radiation  Pattern  of  Surface  Waves  from  Point  Sources 
in  a  Multilayered  Medium,  Bull.  Seismol.  Soc.  America,  54,  377-393. 

McGarr ,  A,  1972,  Comments  on  some  Papers  Concerning  Amplitudes,  Seismic 
Surface  Waves,  J.  Geophys.  Res.,  77,  3823-3826. 

Mendiguren,  J.  ,  1971,  Focal  Mechanism  of  a  Shock  in  the  Middle  of  the  Nazca 
Plate,  J.  Geophys.  Res.,  76,  3861  -3879. 

Salto,  M.  ,  1967,  Excitation  of  Free  Oscillations  and  Surface  Waves  by  a  Point 
Source  in  a  Vertically  Heterogeneous  Earth,  J.  Geophys.  Res.,  72, 
3689-3699. 


Texas  Instruments  Incorporated,  1970a,  Long-Period  Experiment  Software 
Package,  Services  Group,  Dallas,  Texas. 

Texas  Instruments  Incorporated,  1970b,  Array  Evaluation  Software  Package, 
Services  Group,  Dallas,  Texas. 


VI- 1 


Tryggavasson,  E.  ,  1965,  Dissipation  of  Rayleigh  Wave  Energy,  J.  Geophys. 
Res.  ,  70,  1449-1455. 

Tsai,  Y.  B.  ,  and  K.  Aki,  1970a,  Source  Mechanism  of  the  Truckee,  Califor¬ 
nia  Earthquake  of  September  12,  1966,  Bull.  Seismol.  Soc.  America, 
60,  1199-1208. 

Tsai,  Y.  B.  ,  and  K.  Aki,  1970b,  Precise  Focal  Depth  Determination  From 

Amplitude  Spectra  of  Surface  Waves,  J.  Geophys.  Res.,  75,  5729-5743. 

Tsai,  Y.  B.  ,  1972a,  Utility  of  Tsai's  Method  for  Seismic  Discrimination; 

Semiannual  Technical  Report  No.  1,  Texas  Instruments  Incorporated, 
Dallas,  Texas. 

Tsai,  Y.  B.  ,  1972b,  Utility  of  Tsai's  Method  for  Seismic  Discriminatioi  ; 

Semiannual  Technical  Report  No.  2,  Texas  Instruments  Incorporated, 
Dallas,  Texas. 


VI-2 


